Astronomy & Astrophysics manuscript no. prox 


© ESO 2008 


February 5, 2008 





A possible activity cycle in Proxima Centauri 



C. Cincunegui^ ^'*, R. F. Diaz^"^ *, and P. J. D. Mauas 



1,2,* 



,1,3,* 



' Institute de Astronomia y Fi'sica del Espacio, C.C. 67 Sue. 28 (1428) Buenos Aires, Argentina 
2 Fellow of the CONICET 

^ Member of the Carrera del Investigador Cientifico, CONICET 



o 
o 

(N 



o 



> 

in 
m 
o 

o 

Oh 

6 



Received jdate^; accepted [date^ 



ABSTRACT 



Context. Several late-type stars present activity cycles resembling the Solar one. This fact has been observed mostly in stars ranging 
from F to K, i.e., in stars with a radiative core and an outer convective layer. 

Aims. This work aims at studying whether an activity cycle can be detected in the dM5.5e star Proxima Centauri, which is supposed 
to be completely convective. 

Methods. We present periodical medium-resolution echelle observations covering the complete visual range, which were taken at 
the CASLEO Argentinean Observatory. These observations are distributed over 7 years. We discarded the spectra that present flare 
activity, and analyze the remaining activity levels using four different statistical techniques to look for a period of activity. 
Results. We find strong evidence of a cyclic activity, with a period of ~442 days. We also estimate that the Ca II S index varies around 
130% due to activity variations outside of flares. 

Key words. Stars: activity - Stars: chromospheres - Stars: flare - Stars: individual: Proxima Centauri 



1. Introduction 

It is generally accepted that magnetic activity in late-type stars, 
and in particular activity cycles like the one observed in the Sun, 
are the product of the cQ dynamo, which results from the action 
of differential rotation at the tachocline (the interface between 
the convective envelope and the radiative core). Therefore, the 
presence and characteristics of activity cycles are closely related 
to the existence and depth of an outer convection zone. Since this 
depth depends on spectral type - from F stars that have shallow 
convection zones to middle M stars that are totally convective -, 
it is of special interest to study these cycles in stars of different 
spectral types, and in particular in middle-M stars, to determine 
whether there is an onset of cyclic activity. 

To date, activity cycles have been detected i n several late- 
type stars (see, for example, iBaliunas et al.lll985h . mainly mea- 
suring variations in the well known S index, essentially the ratio 
of the flux in the core of the Ca II H and K lines to the contin- 
uum nearby. This index has been defined at the Mount Wilson 
Observatory, where an extensive database of activity has been 
built over the last three decades. However, these observations 
are mainl y concentrated on sta rs ranging from F to K (see, for 
example. IBaliunas et al.|[T995l) . due to the longexposure times 
needed to observe the Ca II lines in the red and faint M stars. 
As a contribution to the subject, we explore the existence of an 
activity cycle in the dM5.5e star Proxima Centauri, using obser- 
vations we made at the CASLEO Observatory, in Argentina and 
that span 7 years. 
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The triple star system a Cen is located at 1.3 pc from 
the Sun. One of its members is Proxima Centauri (a Cen C; 
GJ551), the closest star to Earth, which has a visual magni- 
tude of 11 .01 . Prox Ce n has an effective tempe rature of 3042 K 
dSegransan et alj|2003l) . a ra dius of ~_\J1Rq (lPettersenlll980h . 
and a mass of ~ 0.2 jAlIenl Il976l) . Its rotation period 
has b een estimated several times. From photometric o bserva- 
tions. iBenedict et al.l ( ll998al) found f rot=83.5 days, while lDovl3 
^7^ 



('1987') measu red f rot=5 1 + 12 days from chromospheric activity 
and Guinan & Morgaiil (Il996l) . using lUE observations of h and 
k Mg II lines, found Prot=3 1.5 + 1.5 days. These rotation periods 
are longer than in others M stars (see, for example, Barnes 20031) . 
This longer period is consistent with the older age of Prox Cen. 
If Proxima is co eval with its two companions, its age should be 
4-4.5 10'^ years dPemarque e t al. 1986). 

Prox Cen is also a very active star: it is the first star 
where evidence of a st ellar corona outside flares was found 
jHaisch & Linskvlll980l) . S everal flares have been observed in 
this star, not only in X rays ('Reale et al.' 1988,'2004| ;lGudel et all 
2002, 2004), but also in UV (Haisch et al. 1977), and simultane- 
ously in both (Byrne & McKay 1989; Haisch et al. 1990) . Some 
photometric a nd spectro scopic flares have bee n detected ( i Walkeii 
119811: lPatteni ri994: Ben edict et al]ll998bl) and lSlee etalJ(l2003h 
det ected flares in radio waves. 

iHaisch et al.l (IT990) suggested th at Proxima presents an ac- 
tivity cycle. Using UV observations. iGuinan & MorganI (Il996h 
observed that short-term activity variations have a minimum 
around 1995, and deduced from this fact that the activity cy- 
cle had a minimum at that year, alth ough during that time four 
flares were detected. IBenedict et alj ( Il998a) deduced the exis- 
tence of an activity cycle with a period of about 1 100 days from 
photometric observations . 

This paper is organized as follows: In Sect. |2] we present our 
observations; in Sect. [3] we discuss possible activity indicators. 
We analyze whether our activity measurements are compatible 
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Table 1. Log of the observations. Column 1: a label used in the figures; Col. 2: xJD=JD-2451 000, where JD is the Julian date; 
Col. 3: exposure time in minutes; Col. 4: spectrum clean of cosmic rays ('y') or not ('n') (whether it can be used in the Ca II and 
Ha windows); spectra indicated with '-' do not include the corresponding region, those with are the ones not used in Sect.|4] 
due to the presence of flares. 



label 


xJD 




Ca/Ha 


0399a 


241.85 


25 


-ly 


0399b 


241.87 


25 


-/y 


0399c 


241.89 


15 


-In 


0399d 


242.67 


30 


yh 


0399e 


242.70 


30 


yh 


0300a ★ 


626.82 


30 


yly 


0300b ★ 


626.84 


30 


yly 


0300c ★ 


626.86 


30 


yly 


0300d ★ 


626.88 


30 


v/v 


0300e 


627.82 


45 


yly 


0300f 


627.86 


45 


yly 


0301a 


972.78 


60 


yly 


0301b 


972.83 


60 


yly 


0301c 


974.70 


90 


yly 


030 Id 


974.77 


90 


yly 


0701a 


1 096.46 


90 


yly 


0701b 


1 096.53 


90 


yly 


0302 


1 364.69 


45 


yly 


0602a 


1451.66 


45 


yly 


0602b 


1451.69 


45 


yly 


0802a 


1519.50 


45 


yly 


0802b ★ 


1519.53 


45 


yly 



label 


xJD 


/ 


Ca/Ho' 


0303a ★ 


1715.72 


45 


y/y 


0303b 


1715.75 


45 


yly 


0303c 


1 716.60 


60 


yly 


0303d 


1716.65 


60 


yly 


0303e 


1716.69 


60 


yly 


0303f 


1716.74 


60 


yly 


0303g 


1716.78 


60 


yly 


0303h 


1716.82 


60 


yly 


03031 


1 716.87 


60 


v/v 


0603a 


1 804.59 


45 


n/n 


0603b 


1 804.63 


45 


n/n 


0903a 


1 895.51 


33 


y/y 


0903b 


1 895.53 


33 


yly 


0604a 


2 160.59 


45 


n/y 


0604b 


2 160.63 


45 


n/y 


0604c 


2 161.47 


60 


n/n 


0604d 


2 161.51 


60 


y/y 


0604e 


2 161.55 


60 


y/y 


0604f 


2 161.61 


60 


yly 


0604g 


2 161.65 


60 


yly 


0604h 


2 161.69 


60 


yly 


06041 


2 161.73 


60 


n/n 



label 


xJD 


t 


Ca/Ha 


0904a 


2 274.48 


30 


y/y 


0904b 


2 274.50 


30 


y/y 


0305a 


2448.78 


45 


y/y 


0305b ★ 


2448.81 


45 


y/y 


0605a 


2 523.55 


60 


y/y 


0605b 


2 523.60 


60 


y/y 


0605c 


2 523.64 


60 


y/y 


0605d 


2 523.69 


60 


y/y 


0605e 


2 523.74 


60 


y/y 


0605f 


2 523.78 


60 


y/y 


0605g 


2 543.60 


90 


n/y 


0605h 


2 543.66 


90 


n/y 


0805a 


2 600.51 


80 


y/y 


0805b 


2 600.57 


80 


y/y 


0206a ★ 


2 780.80 


80 


y/y 


0206b 


2 780.86 


80 


y/y 


Total 






49/53 




6200 



Fig. 1. In this figure we show two of the spectra with the best signal-to-noise relations for this star Each of them is the combination 
of seven spectra with an exposure time of one hour The upper spectrum was obtained in March 2003 (xJD 1716) and the one in the 
bottom in June 2004 (xJD 2161). They are arbitrarly displaced for clarity. It can be seen that the spectrum at the top shows a larger 
activity level. 



with cyclic variations in Sect. |4l Finally, we summarize and dis- 
cuss our results in Sect.|5] 



2. The observations 

Our observations were made at the 2.15 m telescope of the 
Complejo Astronomico El Leoncito (CASLEO), which is lo- 
cated at 2552 m above sea level, in the Argentinean Andes. 
We used a REOSC spectrograph to obtain medium-resolution 
echelle spectra covering the wavelength range 3900-7000 A at a 
spectral resolution of K-Aj l^A «26400. We process the observa- 
tions according to the method outlined in Cincunegui & Mauas 
( |2004|) . obtaining flux-calibrated spectra. In Table [T] we show 
the observation logs of Prox Cen. The first column includes 



a label used in the figures, built from the month and year of 
the observation, and a letter to distinguish between different 
spectra from the same observation run; the second column lists 
xJD=JD-2451 000, where JD is the Julian date at the begining 
of the observation and the third column the exposure time (in 
minutes) of each observation. We inspected the regions of inter- 
est of each spectrum to determine if they were clean of cosmic 
rays: in the last column we indicated with 'y' or 'n' whether the 
spectra is suitable or not to be used in the windows coiTespond- 
ing to Ca II and Ha, as explained in the next section. With '-' we 
indicate some spectra that do not include the region of interest 
due to a different instrumental configuration. 

There is a total of 60 individual observations, which have 
been carried out on 24 nights distributed over 7 years. In partic- 



C. Cincunegui, R. F. Diaz, and P. J. D. Mauas: A possible activity cycle in Prox Cen 



3 



ular, on two of the nights we observed Proxima for seven hours, 
obtaining spectra with very high signal-to-noise ratios for this 
faint star, which are shown in Fig. [T] The upper spectrum cor- 
responds to the observations on the night of March 17th, 2003 
(spectra 0303c to 03031), and the lower one was obtained the 
night of June 6th, 2004 (spectra 0604c to 06041). Although no 
flares occurred on either night, as we will show, the overall ac- 
tivity level of the first one is much higher than that of the second 
one, as can be noted by the intensity of several chromospheric 
lines, in particular the Ca II and the Balmer lines, seen in emis- 
sion in this dMe star. 



3. Activity measurements 

The widely-used S index from the Mount Wilson Observatory is 
the most common and best known chromospheric activity indi- 
cator for late-type stars. It is obtained from measurements taken 
with a four-channel spectrometer: two of the channels are cen- 
tered in the cores of the H and K Ca II lines, with a triangular 
profile of 1 .09 A FWHM, and the other two channels acquire the 
flux in the continuum nearby, in two passbands 20 A wide. The 
S index is built as the ratio between the chromos pheric H and 
K flux es and the photospheric continuum fluxes (IBaliunas et al.l 
119981) . 

To reproduce Mount Wilson's S index we integrated the 
fluxes in our spectra around the K and H Ca II lines, weighting 
them with a triangular profile resembling the instrumental re - 
sponse of the Mount Wilson spectrometer (iVaughan et al.lll978l) . 
We normalized these fluxes by the average flux in the same two 




Fig. 2. Comparison of the Ca II index, S , and the Ho- index. A, 
for the 47 spectra in which both spectral features have been mea- 
sured. The errors are assumed to be 12% in A and 17% in S . 
The least-squares fit has a slope of 0.19 ± 0.03 and a correla- 
tion coefficient of -0.94. The fit significance is 60%. The dashed 
lines represent the points that deviate +3cr from the fit. In the 
small figure the confidence regions of 39.3% (full line) and 90% 
(dashed line) are drawn. 



20 A continuum passbands used at Mount Wilson, centered at 
3901 and 4001 A. In a previous work ( Cincunegui & MauasI 
2002), we checked the accuracy of our m easurements observ - 
ing a set of 18 non-variable stars given bv lHenrvet all ( 119961) . 
and using them to intercalibrate our measurements with Mount 
Wilson's S. We have found a very good correlation between 
both sets. However, Proxima Centauri is a very red a nd faint star 
{B - y=1.807 and V^l l.Ol. lPerrvman et al.lll997l) . Therefore, 
the Ca II lines are very hard to observe with a good signal-to- 
noise ratio, requiring very long integration times with a 2m class 
telescope. 

To overcome this problem, it is worth noting that in Proxima 
Centauri all the other Balmer lines, and in particular Ha, are ob- 
served in emission, as is the case for the most active late-type 
stars, the dKe and dMe stars. Therefore, we constructed an in- 
dex A from the H„ profile, defined as the ratio between the av- 
erage flux in a 1.5 A square passband centered in the line and 
the average flux in a 20 A nearby continuum window, centered 
at 6605 A, in a similar way to the definition used for S . 

To check whether A can be used as an activity index for 
this red star, we explored the relation between A and our S in- 
dex, which is shown in Fig. |2] We estim ated the errors in the 
flux-calibration method employed at 10% (Cincunegui & MauasI 
i2004). Therefore, 10% should be considered as a lower limit for 
the error in each index. Using a non-linear minimization to 
fit a straight line considering the errors in both coordinates, and 
varying the percentual errors for each index until finding a rea- 
sonable significance of 60%, we found a very good correlation 
between both indexes, with a correlation coefficient of 0.94 and 
a slope of 0.19 for errors of 17% in S and 12% in A. The dif- 
ference between the errors is natural since the Ca II lines have a 
much lower signal-to-noise ratio. In the insert to Fig.|2] we show 
the confidence regions for the parameters of the fit. In spite of 
the much lower sensitivity of H,,, as indicated from the small 
slope, we chose the A index to study activity in this star because 
of the much better signal-to-noise ratio. 

Since Proxima is a flare star, it is expected to have frequent 
flares. Therefore, before exploring the existence of an activity 
cycle in Proxima Centauri, we have to determine whether any 
flare occurred during our observations, since flare activity can 
mask the variations of activity due to the cycle, and it has to be 
filtered out before the analysis. In Figs.|3]and|4]we show the line 
profiles of the Ca II H and K and Ha lines. Visual inspection 
of the temporal series of observations reveals the possible oc- 
currence of several flares. The first one corresponds to the four 
observations 0300a to 0300d: the apparent flare goes from max- 
imum to minimum, with the two central observations with sim- 
ilar levels of activity. Also, in the spectrum 0802b, all the lines 
are much stronger than in the previous one 0802a. Similar varia- 
tions can be appreciated between observations 0303a and 0303b, 
0305b and 0305a, and 0206a and 0206b, although in these cases 
the lines are not so intense. Using these considerations, we sep- 
arated the flaring spectra from the "normal" ones. We excluded 
these flaring spectra - which are indicated in Table[T]with a '★' - 
from the rest of the analysis. 



4. Behavior of l-ia with time 

For the non-flaring spectra, we calculated a nightly average of 
the A index, shown in Table |2] which we examine in this section 
to determine the possible presence of an activity cycle. Using 
this data, we first calculated the Lo mb-Scargle periodogram 
(IScargldll982L lH orne & Baliunas|[T986i) using the algorithm in- 



4 



C. Cincunegui, R. F. Diaz, and P. J. D. Mauas: A possible activity cycle in Prox Cen 



Table 2. Nightly averaged A index. Column 1 : the label used 
in Fig. |7j Col. 2: the xJD at the beginning of the observations; 
Col. 3: the quantity of spectra averaged; Col. 4; the total expo- 
sure time (in minutes); Col. 5: the average index. 



Labe: 




xJD 


A' 




A 


0399 




241.85 


2 


50 


2.89 


0300 




627.82 


2 


90 


2.50 


0301. 


1 


972.78 


2 


120 


2.14 


0301. 


2 


974.70 


2 


180 


2.00 


0701 




1 096.46 


2 


180 


2.84 


0302 




1 364.69 


1 


45 


2.16 


0602 




1451.66 


2 


90 


2.13 


0802 




1519.50 


1 


45 


2.84 


0303. 


1 


1715.75 


1 


45 


2.28 


0303. 


2 


1716.60 


7 


420 


2.63 


0903 




1 895.51 


2 


66 


2.46 


0604. 


1 


2 160.59 


2 


90 


2.35 


0604. 


2 


2161.51 


5 


300 


2.00 


0904 




2 274.48 


2 


60 


2.48 


0305 




2448.78 


1 


45 


2.61 


0605. 


1 


2 523.55 


6 


360 


2.46 


0605. 


2 


2 543.60 


2 


180 


2.54 


0805 




2 600.51 


2 


160 


2.66 


0206 




2 780.86 


1 


80 


2.42 



eluded in iPress et al] (|1992|) . This periodogram is a method to 
estimate the power spectrum (in the frequency domain), when 
observation times are unevenly spaced, and normalizes the spec- 
trum in such a way that it is possible to measure the signifi- 
cance of the peaks. The amplitude of the periodogram at each 
frequency point is identical to the equation that would be ob- 
tained estimating the harmonic content of a data set, at a given 
frequen cy, by linear least- squares fitting to an harmonic function 
of time dPress et al.lll992l) . We show this periodogram in Fig.|5] 
It can be seen in the figure that there is a distinct peak at 442 
days, with a significance of 65%. In other words, a period of 442 
days is detected, with a False Alarm Probability (PAP) of 35%. 
It is worth noting that this FAP is very good, since the activity 
cycles are neither harmonic nor exactly periodic. Other less sig- 
nificant peaks are present, the most important one corresponding 
to 166 days with a FAP of 80%. This peak could correspond to a 
subharmonic of a year, which can appear due to the fact that our 
observations are carried out with a periodicity of approximately 
one year. 

We also applied several techniques in the time domain to 
check the existence of this period. Basically, all of them start 
with the construction of the light curve (that is, the observed 
data as a function of phase for different trial periods). Then, this 
light curve is partitioned in several bins, whose quantity depends 
on the amount of data available. In the three cases, we chose a 
partition of two bins, due to the small number of data points. 

Th e first of these techniques is the one outlined in lJurkevichI 
( [1971 . which assumes that the real period is the one that mini- 
mize the dispersion, measured from the avera ge, of t he combined 
bins of the light curve. The second one jMarraco & Muzzid 
Il980l) is an improvement of the first, which consists in measur- 
ing the dispersion relative to the leas t-squares straight line for 
each bin. Finally, ICincotta et al.l ( 1199 5) instead of minimizing 
the dispersion, calculated the so-called Shannon entropy, which 
is a measure of the order of the data. Again, the real period will 
minimize this quantity. 



In Fig. |6]we show the dispersion of the light curve - or the 
Shannon entropy - as a function of the trial period for each 
method. In all of them we find periods compatible within 10% 
with the one of ~442 days (420, 492, and 412 days). Also present 
are possible harmonics of a year (198, 359, and 1090 days). One 
should keep in mind that solar activity is not exactly periodic, 
since each cycle has a different length and intensity. Therefore, 
we believe that the agreement found between diff'erent methods 
is quite remarkable. 

We also fit the data, by least-squares, with a harmonic func- 
tion with a period of 442 days, and we obtained the curve shown 
in Fig.lTJ;. Also the fit is remarkably good here, again consider- 
ing that the cycle is not exactly harmonic. Finally, in Fig.lTJ? we 
draw the associated light curve. 

As measured from minimum to maximum, the amplitude of 
the period is around 25% in A. If we translate this amplitude to 
the S index, using the slope of Fig.|2] we find a variation of about 
130% in a c ycle, to be cornpared with a solar variation for S of 
only ~30% dBaliunas et al.|[T995h . This much larger variability 
is consistent with the fact that Proxima Centauri is an extremely 
active star Since there are no sy stematic studies of ac tivity varia- 
tions for M stars so far (e.g., see lBaliunas et al.ll995l) . we cannot 
compare with any other similar star 

Finally, we checked the accuracy of the result of Fig. fus- 
ing the sunspot numbers taken from the National Geophysical 
Data Center. To do so, we took a sample of the solar data with 
the same phase intervals that we have in our data for Proxima, 
we added Gaussian noise with errors of 10% at each point, and 
we computed the Lomb-Scargle periodogram. We repeated this 
procedure 1000 times with random starting dates. 

For each periodogram, we considered the period with maxi- 
mum significance (or minimum FAP) as the detected period . In 
Fig. |8]we plot the histogram of these detected periods. The full 
line shows only the periods with FAP < 0.35. We see that the 
correct period (10 to 12 years) was detected in 62% of the peri- 
odograms and in 44% of the cases with a FAP < 0.35. Only in 
6% of the cases was a "false" period detected with a FAP < 0.35. 

period [day] 
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Fig. 5. Lomb-Scargle periodogram of the data of Table |2l The 
False Alarm Probability levels of 50 and 80% are shown. 
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) 


0.9 




0.85 




rfal98 days 



-r^420 days 
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Fig. 6. The results of the three different methods employed 
that involve light-curve dispersion measurements: in the higher 
panel, the relative Jurkevich (1971) dispersion; in the mid- 
dle panel, the iMarraco & Muzzio ( 1980) dispersion; and in the 
lower panel, the Shannon entropy CCincotta et al.„1995i) . for the 
data ofTableE] 



Therefore, if this star has a cyclic behavior similar to the 
solar one, the probability of detecting it in our observations is 
P ~ 60% and the probability of detecting it with a FAP < 0.35 
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Fig. 7. fl) A index vs. time, assuming errors of 5% for each point. 
The full curve is the least-squares fit to the data of Table |2] after 
establishing a harmonic function of period 442 days. The dotted 
curves represent the +2cr deviations, b) Light curve of the same 
data. 
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Fig. 8. Histogram of the detected periods (the periods for which 
we found the maxima Pls) for a sample of 1000 different sets 
taken from the Sun. The full line shows only the periods with 
FAP < 0.35 and the dotted line shows all the cases. 



is P ~45%. On the other hand, if a period is detected with a 
FAP < 0.35, its value is correct with a probability of P~90%. 



5. Discussion 

In this paper we found strong evidence of a cyclic activity in the 
dMe star Proxima Centaury, with a period of ~442 days +10%. 
Similar values for the period were found from the Lomb-Scargle 
periodogram and using three different techniques in the time do- 
main. We were also able to determine that the activity variations 
outside of flares amount to 130% in S, four times larger than for 
the Sun. 

As we mentioned before, the possible existence of an activ- 
ity cycle in th is star is mentioned in three previous papers. First, 
iHaisch et al.l (ll990) found some evidence su^ ^estin^ a stellar cy- 
cle, al though they did not estimate its period. iGuinan & MorganI 
(Il996h found some manifestation that in 1995 Proxima was near 
a minimum of activity, again without a n estimate of its period. In 
a more detailed study, .Benedict et al.l ( il998a) found some hints 
of an activity cycle of about 1100 days, compatible with the 
minimum activity around 1995, although they had some gaps in 
temporal coverage. The middle point of what they consider min- 
imum activity corresponds to JD 2 449 750. If we extrapolate our 
period back in time, taking into account that in JD 2 45 1 500 we 
have minimum activity, we find that in J D 2 449 732 we a l so have 
minimum activity, very close to what [Benedict et al.l ( Il998ah 
found. However, the data used in these two last works were also 
used to stu dy the rotational per iod of this star, giving very differ- 
ent results: iGuinan & Morganl (11996.) found f iot=31.5+1.5 days 
and [Benedict et al. (Il998al) det ermined Prot=83.5 days. On the 
other hand, Kiirste r et al. f ll999l) did not find any evidence for the 
rotation periods found in these two works, based on very precise 
radial velocities of Proxima Centauri, which they analyze using 
both the Scargle periodogram and a sine-fitting routine, which 
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ininimize taking into account data errors. Therefore, we be- 
lieve that our value of ~442 days for the cycle is more reliable, 
in particular because all the methods we employed give similar 
values for the period, and because we estimated, using solar data, 
that if a cycle is present, the period we found is the correct one 
with a probability P ~ 90%. 

Concerning the flare frequency, IWalkeij (Il98ll) predicts a 
flare of intensity greater than 5 10^^ in the U band once ev- 
ery 0.8 hr For the most intense flares, the predicted frequency 
falls to one every 3 1 hr In our case, we observed Proxima for 
~54 hr, detecting two very strong flares (0300a to d and 0802b) 
and three weaker ones (0303a, 0305b, and 0206a) in this period. 
Therefore, we found a frequency of one flare every 10 hr 

It is generally accepted that magnetic activity results from 
the generation of a large-scale toroidal field by the action of dif- 
ferential rotation on a poloidal field, at the interf ace between 
the convective envelo pe and the radiative core ( PaAe3[T97l 
ISpiegel&Zahnll992l) . This is the o-O dynamo, which predicts a 
strong correlation between activity and rotation, which is in fact 
observed for stars from F to early M. This theory also explains 
the fact that young rapid rotators - which have larger differential 
rotation - have higher activity levels than old slower rotators, 
with less differential rotation and therefore a less important dy- 
namo eff'ect. 

On the other hand, stars ranging from M3 to M9 are thought 
to be fully convective, and the refore th ey do not support an aO. 
dynamo. It has been argued dBarne sl l2003h that fully convec- 
tive stars do not possess large-scale dynamos, and that their ro- 
tational morphology are dominated instead by small-scale tur- 
bulent fields. Nevertheless, there is plenty of observational ev- 
idence that slow late-type rotators like dMe stars have strong 
magnetic fields, with filling factors larger than for earlier stars 
and enhanced activity (^ Mochnacki & Zirin|[T980l;lHawlevll 19891: 
iMauas & Fal chi 1996)^ 

Recently, Cha brier & Kiikeij (2006) showed that these cool 
objects can support large-scale magnetic fields by a pure dy- 
namo process. In this dynamo, helicity is generated by the ac- 
tion of the Coriolis force on the convective motions in a rotat- 
ing, stratified fluid, and produces an amplification of the mean 
magnetic field. Moreover, these fields can produce high levels of 
activity like the ones observed in M stars. This dynamo does 
not predict a cyclic activity. However, our observations suggest 
that this cool star does have a clear period. 

The activity cycles measured in earlier stars (F to K) are all 
longer than 2.5 years ( Baliunas et al.lll995l) . to be compared to 
the one we found for Proxima Centaury, of only 1.2 years. On 
the other hand, the variation of the activity levels is much larger 
for this star than for earlier-type ones. As we explained, the dy- 
namo models show that the generation of activity in M stars is 
a very different phenomenon than in earlier stars, i.e., stars with 
a radiative core and an outer convection zone. It is desirable to 
confirm the existence of a cycle in Prox Cen and its period, and 
to explore other M stars near the limit where stars become fully 
convective, to constrain the dynamo at work in these stars. 
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Fig. 3. Ca II region for each 
observation of Table [T] In each 
window the horizontal axis 
goes from 3890 to 4012 A, and 
the vertical one from -0.05 to 
1.2 10^'^ erg cm-- s"' A"'. 
Some spectra (0603a and b; 
0604a, b, c, and i; and 0605g 
and h) are deficient and are not 
considered in the subsequent 
analysis. The first three spectra 
were acquired with a different 
instrumental configuration and 
they do not include this region. 
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Fig. 4. Ho- region for each ob- 
servation of Table [T] In each 
window the horizontal axis 
goes from 6540 to 6620 A, and 
the vertical one from 0.15 to 
2.25 10-'2 erg cm-^ s"' A"'. 
Some spectra (0399c, 0603a 
and b, and 0604c and i) are de- 
ficient and are not considered 
in the subsequent analysis. The 
fourth and fifth spectra were 
aqcuired with a different instru- 
mental configuration and they 
do not include this region. 
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